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Abstract 

In a previous article a relationship was established between the 
linearized metrics of General Relativity associated with geodesies and 
the Dirac Equation of quantum mechanics. In this paper the exten- 
sion of that result to arbitrary curves is investigated. A generalized 
Dirac equation is derived and shown to be related to the Lie deriva- 
tive of the momentum along the curve. In addition, the equations 
of motion are derived from the Hamilton-Jacobi equation associated 
with the metric and the wave equation associated with the Hamilto- 
nian is then shown not to commute with the Dirac operator. Finally, 
the Maxwell-Boltzmann distribution is shown to be a consequence of 
geodesic motion. 

KEY WORDS: non-geodesic motion, Dirac equation, equations of 
motion, Maxwell-Boltzmann distribution. 



1 Introduction 

In a recent paper Frank Tipler [10] gives a derivation of the Schrodinger 
equation using the Hamilton-Jacobi principle of action. In doing so, he 
is able to transpose Erwin Schrodinger's "purely formal procedure" [1] of 
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replacing in the Hamilton- Jacobi equation with ±7^^ with a more 
direct derivation of the wave equation. Likewise, in a paper by Marie- 
Noelle Celerier and Laurent Nottale [4] in 2003, the Dirac equation is de- 
rived by taking the square root of the Klein-Gordan equation using the 
"bi-quaternionic action" associated with geodesic motion and the Hamilton- 
Jacobi principle of action. Other approaches not directly based on the use 
of the action principle have also been used. For example, Ng and Dam used 
a geometrical derivation of the Dirac equation by exploiting "rotational in- 
variance" and "the explicit use of the spin-| property of ip" [5], while Martin 
Rivas quantized a Poincare invariant Hamiltonian in which the spin angu- 
lar momentum of the particle is constant with respect to the center of mass 
observer to derive the Dirac equation. [9] 

This paper contains parallels with the approach of Tipler, and Celerier 
and Nottale in that it makes use of the Hamilton Jacobi equation but it 
also differs in that it emphasizes the fundamental role of the metric in 
enabling us to derive two equations associated with non-geodesic motion: 
one a generalized Dirac equation, which captures the kinematics, and the 
other an equation of motion describing the dynamics. 

Throughout the paper, (Ai,g) will denote a connected four dimensional 
Hausdorff manifold, with metric g of signature -2. At every point p on the 
space-time manifold Ai we erect a local tetrad eo(p), ei(p), ^(p), e-sip) such 
that a point x has coordinates x = (x°, x 1 , x 2 , x 3 ) = x a e a in this tetrad 
coordinate system, while the spinor \1/ can be written as ^ = ^V^p), where 

represent the coordinates of the spinor with respect to the tetrad at p. 
This is permissable since a spinor is also an element of a vector space. In 
particular, when each ip % is equal for each % then we can write \I/ = ipt; where 
the spinor £ = (e$) and ||\&|| = 1. When this occurs, we say that ip is a 
scalar field. Such a scalar field will occur when \& is parallel transported 
along a geodesic. For example, in the case of the Dirac equation, the plane 
wave solution can be written as 

V(x) = 

where ip = e ~ ipx l h and u + (p) 
solutions for each momentum p 



- e- l ^ h u + (p), (1) 

= f Ul I and allows for two independent 
. Also in this case, the spinor u(p) remains 
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constant along a geodesic. This reflects Rivas result in that the derivation 
of the Dirac equation presupposes the motion of a particle "in a plane 
orthogonal to S, which is constant in this frame," [9] where S refers to spin- 
angular momentum. In terms of quantum mechanics each constant can be 
identified with a quantum number. Also at each point p we can establish 
a tangent vector space T p (Ai), with basis {do, d\, 82, 83} and a dual 1-form 
space, denoted by T* with basis {dxo, dx\, dx2, dx^} at p, defined by 

dx"d v = d v x» = 5%. (2) 

We refer to the basis {dx°, dx 1 , dx 2 , dx 3 } as the basis of one forms dual to 
the basis {do, d±, <9 2 , d 3 } of vectors at p. 

With notation clarified, we note that in a previous article [8] , the Dirac 
equation was derived as a dual of a linearized metric tensor for geodesies. 
The linkage was accomplished in a natural way by associating a gener- 
alized Dirac equation with those operators which are duals of differential 
one-forms, obtained by linearizing the metric tensors of General Relativity 
(expressed locally as a Minkowski metric). Specifically, in that paper, if 

ds 2 = g iW dx iX dx v = r] ab dx a dx b (3) 

where a and b refer to local tetrad coordinates and i] to a rigid Minkowski 
metric tensor of signature -2, it was shown that a matrix ds = , ~f a dx a with 
eigenvalue ds could be associated with the metric tensor ds <8> ds (which 
from now on we will denote by ds 2 ) such that {~i a ,lb} = ^Vab- 

In addition, by noting that the ds matrix is the dual of the expression 
d s = ^ a -^ this enabled us to associate a generalized Dirac equation 

OX a OS 

with the motion of a particle along a geodesic, with a fixed spinor £. 

This paper is an extension of that previous work. First, we re-investigate 
this result by showing that the existence of a Hamilton- Jacobi function asso- 
ciated with a family of space filling curves is equivalent to defining the Dirac 
equation at almost every point on the manifold. Indeed, this association 
also suggests a deeper understanding of the term "wave-particle duality." 
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Secondly, we show that the existence of such functions also determines the 
existence of both quantum and classical ideal gases defined in terms of the 
Bose-Einstein and Boltzmann statistics respectively. Thirdly, we introduce 
accelerations and then show that the integrability of the Hamilton- Jacobi 
function determines a coherent set of natural motions and also implies the 
non-commutativity of the Hamiltonian and action operators associated with 
non-geodesic motion. 



2 Generalized Dirac Equation 



The generalized Dirac equation defined above relies on the definition of the 
four-momentum in special relativity and upon the fact that d s ifj and ds are 
parallel along geodesies, and consequently by the chain rule their product 
is ^jjds. In contrast, when accelerations are introduced we will find that in 
general 



ds sy 2\ds , sV J 2 

dip ds d^> 
ds ds ds 



ds 



(5) 
(6) 



In this regard, note that equation (5) represents an operator equation and 
the four 7 matrices used to define ds = ^ a dx a and d s ip = l a -§^- are a 
spinor representation of four unit vectors constituting a local Minkowski 
tetrad at any point on the curve. Equation (6), on the other hand, is 
written in the usual vector notation. The two equations can be identified 
by noting that the anti-commutator and commutator relationships in (5) 
associated with the spinor tetrad, define a dot product and a cross product 
respectively, which can can also be expressed in standard vector notation 
as in equation (6). Direct calculation of the dot product of ds and its dual 

d s ip or equivalently of dx a and J^. gives ^ds. 

This latter term can also be directly related to the Hamilton- Jacobi 
characteristic function, which in turn can be associated with a coherent set 
of natural motions [3]. Indeed, to remove any ambiguity, we begin with the 
following definitions: 
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Definition 1 A function W = J a ^ (p^f — H-^j dX is called a Hamilton- 
Jacobi function if the integral is path independent for all curves in {<r(A) G 
(M,g)} and ^ = -H(x u x 2 , x 3 , t, g, g, g) de/ined respect to a 
local tetrad . 

Equivalently, we can say that dW = pdx — Hdt is an exact differential. 

We begin by focusing on the kinematics associated with the function W 
and by showing how it can be related to the Dirac equation. This leads to 
the following lemmas and corollaries: 



Lemma 1 Let ip(W(x a )) be a differentiate function of W and let ip' = 
W = — f rj ab p a dxb = J pdx — Hdt is the Hamilton- J acobi function, iff 
ip{W) = f p* dx — H*dt is also a Hamilton- J acobi function. 

Proof:(=^) First note p* a = ip'p a , from the differentiability of *p(W), and 
p = H . Since W is a Hamilton- Jacobi function: 

djAAW) dW 

dt ~ ^~dT 

= -ip'H(x 1 ,x 2 ,x 3 ,t,p 1 ,p 2 ,p 3 ) 

= -^'H{x u x 2 , x 3 , t, pl/ip', p* 2 / ip', P* 3 /ip') 

= -H*(x 1 ,x 2 ,x 3 ,t,p* 1 ,p 2 ,p* 3 ). 



Therefore, ip(W) is a Hamilton- Jacobi function. 



(<^=) Conversely if ip(W) is a Hamilton- Jacobi function then 

dW dip(W) 

~df ~ dt 

= -H*(x 1 ,x 2 ,x 3 ,t,p* 1 ,p* 2 ,p* 3 )/ip' 

= -H(xi, x 2 ,x 3 , t, p\l ip', p\l ip', p* 3 /ip') 

= -H(x 1 ,x 2 ,x 3 ,t,p 1 ,p 2 ,p 3 ). 

Therefore, W is a Hamilton- Jacobi function. 
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In addition, from equations (5) and (6) we obtain the following Lemma 
which gives a necessary and sufficient condition for the existence of such 
Hamilton- Jacobi functions. 

Lemma 2 Let ^>(W(x, t)) be a differentiable function and {c(A)} a family 
of curves on the manifold with unit tangent vectors with respect to a local 
tetrad then ds.d s ip{W) is an exact differential iffip(W) is a Hamilton-. Jacobi 
function such that p* a = g^f = where ip(W) = ip\\(W) + *p±(W). 

Proof: From equations (5), if ^.d s ip(W) is an exact differential then 

[f s , ^ff 1 } is also an exact differential. This means that if>(W) = i/>\\(W) + 

ip±(W) where [J, ^p-} = and ip ± (W) = c Q t + c lXl + c 2 x 2 + c 3 x 3 , with 
c , ci, C2, C3 being constants. 

Also, , 9 ^qY^ \ = means ^-is parallel (for any parameter A including 
the curve length s) to ^g^, and 3g(X) such that p\ a = g(X)^ = ^0^- 

Also given ^.d s ip\\(W) is an exact differential and denoting x = t, gives 

'kdtb(W) d ^\\( w ) dxl ] d^ {l (W)dx 2 | <9^|(HQcfo 3 | dj>\\(W) dt 
ds s dx 1 ds dx 2 ds dx 3 ds dt ds 

#j|g) 
ds 

On substituting g(s)^ = ^J, 1 ^ and noting that ^-^ff = 1 gives 
g ( s ) = ?Mp.. It follows that X 

2 



(T) 

Therefore ip\\(W) is a Hamilton- Jacobi function, as also is ip(W) = ip\\(W) + 
i/>±(W) with p* a =p* a + Ca . 
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Conversely, given f a = f ^+c a = «^+^ = then [ds, ^|| 

and since ip(W) is a Hamilton- Jacobi function, it follows from Equation 
(5) and the definition of ip± that 



ds.d,ip\\(W) = dil>\\(W) and 



ds djjjyV) 
ds' ds 



are both integrable. Therefore, ds. d s ijj(W) is an exact differential. The result 
follows. 

Corollary 1 If ip(W) = W = —m s + k is a Hamilton- Jacobi function for 
each fixed k then g(s) = m is called the rest mass associated with linear 
motion in the plane W(s, k), and 

H 2 =pl+pl+pl + m 2 . 
Proof: Follows from Equation (7). 

We now show that for the family of curves defined in Lemma 2, we can 
always construct a Hamilton- Jacobi function locally at each point on the 
curve. 

Corollary 2 Let ds 2 = dx a dx a define a metric tensor locally along a piece- 
wise smooth congruence of curves, {o~(\) = x a (X)}, which fill coordinate 
space with only one curve passing through each point, then a Hamilton- 
Jacobi function can be constructed such that 



f ds i f i dt . , _ dx u d\ 

W — \ m—ds = / pax — m—dt, where p = m——, m = m —. 
Ja dX J dX dX ds 

(8) 

Note ^ = m is called the rest mass. 

Proof: At every space-time point in a gravitational field, by the Prin- 
ciple of Equivalence, an inertial frame can be constructed. Now, define 
p a = m^fr, tu = m(s)4^, where m(s) = m is a constant. Then dW = 

,2 



m 



(if) ^ = m ods is clearly an exact differential, such that 

W = m s + k = m \jx\ + x\ + x\ — t 2 + k, k constant 



in a neighborhood of any point. Therefore d s W = 7 a |^ exists almost ev- 
erywhere and dW = ds.d s ip{W). It follows from Lemma 2 that W is a 
Hamilton- Jacobi function. 

Remark: The corollary has shown the existence of W at every point on the 
manifold by showing that W can be constructed in terms of a local tetrad 
as a linear function of the arc length, s. In effect, this means that we are 
defining the Hamilton- Jacobi function along a geodesic passing through the 
point x . In practice, it is sufficient to integrate |^ = rn ^£- to obtain W, 
as in the example below. Moreover, although p Q = p(x ) is a constant along 
the geodesic passing through x a , this does not mean that p a = |^ is con- 
stant along an arbitrary curve. Indeed, this latter condition only follows, if 
motion is along a geodesic. 

Example: As an example of the corollary, consider a motion of a particle of 
rest mass m fired into the air without resistance. The Equation of Motion 
for such a projectile in Minkowski space are given by 

F = -m gj, (9) 

or equivalently, denoting ^ = x, 

m x = m y = -m g. (10) 

Solving the equations gives 

x = x + u x s and y = y + u y s - ^gs 2 . (11) 

However, in terms of the corollary above, we begin not with equations 
of motion (9) but with the line increments associated with the kinematics 
described by equation (11). Specifically, 

—ds 2 = dx 2 + dy 2 — dt 2 
iff — m ds = m xdx + m ydy — m idt 
iff — m ds = m u x dx + m {u y — gs)dy — m idt. 
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If we now let dW = —m ds along a geodesic and require that W be an 
integrable function passing through (xo,yo,to), we obtain W = m u x x + 
m>o{ u y ~ 9 s )y ~ rn c 2 it + w Q . Indeed, on taking partial derivatives, we find 

dW dW (dW\ 

— = m u x =p x , — = m (u y - gs) = p y \—-\=-m t 

and 

(9W\ 2 2,2,2,2 rj2 

[~Qfj = m o + Pi + P2 + Ps = H ■ 
Therefore, H is a Hamilton- Jacobi function. 

In effect, Lemma 2 establishes a relationship between Hamilton- Jacobi 
functions and the commutator relationship, [ds, d s ifj(W)} applied to Equa- 
tion (5) (or equivalently Equation (6)). We now use the same commutator 
relationship to establish another important property relating the dual op- 
erator d s and the metric operator ds associated with the increments along 
a curve. Indeed, Equation (5) could be described as the most general form 
of a "wave-equation" associated with a curve in space-time. However, we 
put the expression "wave-equation" in quotes to emphasize that ip{W) is 
not necessarily a wave-function of quantum mechanics. For the moment, ip 
they can be any C\ function defined on the manifold. In a previous paper, 
we have noted that ip can only be interpreted as a quantum wave function 
when further restrictions are imposed on the function space, such as requir- 
ing that it be an L 2 function. Intuitively, we could think of dip(W(s)) as 
a wave associated with the vibration of a curve a(s) in space-time, whose 
tangent is ds, with respect to a local tetrad coordinate system. This leads 
to the following lemma: 

Lemma 3 If tp(W) is a Hamilton-, Jacobi function such that [d s W,ds] = 
then there exists a simultaneous eigenfunction £ such that 



(d s ^(p) = d^(p), (12) 

nnnnnn o<n i 

ds 



where dip s = |p ; which in the case of geodesic motion reduces to 



d s ty = — , where V = . (13) 
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Remark: (d s ip)^ = d s ^ in general, since x is independent of p in phase 
space. However, ^- ^ 4>'(p)^ unless motion is along a geodesic. 

Proof: First note that [d s W, ds] = implies [d s ip,ds\ = [ip'(W)d s W,ds] = 
0. Therefore, there exists simultaneous eigenvectors £ = such that 
ds£ = ds£ and (d s ip)£ = 7X£(p) = l a Pa^'£(p) = mcip'(p)£(p) = (d s i/j)£(p). 
Also, is constant along a geodesic and therefore 

dV 

d s ^ = ——, where $ = 
ds 

The result follows. 



Corollary 3 Let ds 2 = dx a dx a define a metric tensor on a manifold and 
p a = m(s)^f- along a curve. Then there exists a Hamilton- Jacobi function 
ip{W), and a vector ^{p) such that 

Proof: By Cor. 2 it is sufficient to define p a = and solve for W to 
obtain a Hamilton- Jacobi equation. It now follows that 

dW 

[d s W, ds] = [ 7 a ^, l a dx a ] = [ 7 V, l a dx a ] = 0. (14) 

But from Lemma (3) if [d s W, ds] = 0, there exists a ip — ip(W) such that 

(4V0£(p) = (d s ^(p). 

The result follows. 



Remark: We refer to Equation (12) as a generalized Dirac equation as- 
sociated with a curve, and Equation(13) as a generalized Dirac equation 
associated with geodesies. It reduces to the usual form of the Dirac equa- 
tion if we let ip = Ac kW , where A is an arbitrary constant and k — I : 

Corollary 4 Let ip = Ae KW where W = —ms + k as defined in Cor. 1 and 



k — 1 then 

n 



i 

7 a — = --m*. (15) 

dx a h v ' 
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Proof: Direct substitution for W in Equation (13) reduces to the conven- 
tional Dirac equation 

7 - — = —-m^. 

ox a n 

The above equation can be rewritten in the conventional form, if we multiply 

across by — i/ry°, define a = 7 , a a = 7°7 a and let E = iTiH* = iTi^- to 
get 



/ d d d \ 

-1T1 an- h «2t^ h a: 3 —— + O! m 

\ 0x1 0x2 0x3 J 



vj> = E^I where * = 

(16) 



Corollary 5 Consider motion in a plane with Hamilton- Jacobi equation 
W = -m Q s + cix 1 + c 2 x 2 + c 3 x 3 + c x° + d = W\\ + W±. Let S a = ^t = c a 
(called the spin), da = ^ a dx a + ^ a (c a /m )ds, d a W = 7 a g^ = >y a (Pa + c a ), 
then there exists an eigenfunction £ such that 

{d^)Op) = (^)e(p). (17) 

Proof: By construction [7 a |^,<i<7] = therefore by Lemma 3 

{d^)Op) = (^)e(p). 

Remark: (1) This corollary confirms Martin Rivas' observation that the 
Dirac equation presupposes the motion of a particle "in a plane orthogonal 
to S, which is constant in this frame," [9] where S is the spin. 

(2) The duality relating line increments and the Dirac equation brings to the 
foreground philosophical issues regarding the difference between quantum 
and classical mechanics, and in particular how the two might be related. A 
more detailed discussion of this point will be found in the next section and 
also in [8]. 

Next, we prove a theorem that relates the Lie derivative and the Dirac 
equation for a particle in a closed system. 
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Theorem 1 Consider a family of curves {o~(s)} parameterized by the curve 
length s with unit tangent vectors u a = ^ associated with the increments 
ds 2 = dx a dx a then the Lie derivative C u (p) = iff there exists a Hamilton- 
Jacobi function ip{W) such that {^ a ^p-, ds] = and 

Proof:(=^) Since C u (p) = there exists a coordinate system (in fact the rest 
frame) [2] such that u a = 5^ and p a Q = 0. Therefore, in a general Lorentz 

frame p a = m u a + m v a where m is constant and = 0. Denote 

m u a by pji. Then pjj = ^ defines a Hamilton- Jacobi function by Cor. 2 

and consequently for tp = ip(W), [d s ip,ds] = 0. It follows from Lemma 3 
that 

(d s mp\\) = (d s mpw)- 

(<=) Recall ^ = 4(W). Therefore J£ = ^'jg. Now [T^Js] = implies 

7 a f| = ™(s)f s for some function m(s). But (d a ip)£(p\\) = (d s ip)^(pn) 
implies m(s) = while Cor. (2) implies that m(s) = m . Now let 
P\\ = f^" = m °^T' -D e fi ne P a = P\\ + m oV a , where u a is a Killing vector for 
v a (i.e.£ u (v) = 0). It now follows that 

C u {p a ) = C u {p\\) + C u {mv a ) 
= (pl)., b u b -u%p\ 
= mu a b u b - mu a b u b 



The theorem has been proven. 



Remark: The theorem states that the Dirac equation associated with a 
particle exists and is defined locally if and only if the wave function is Lie 
transported along the curve whose action is W. 

By way of concluding this section, we make some final observations: 

• In general for any Hamilton- Jacobi function ip(W) it is possible to 
define d s ip = d s ip\\ +d s ip± such that [d s ip\\, ds] = and {d s ip\\, ds} = 0. 
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d s ip\\ is the projected cosine along ds and satisfies 

dW 

= (d s ii))\\Z(p), with Pa = — , 

where ds£(p) = ds£(p). Also, d s ip± defines the spin along ds. 

The Hamilton- Jacobi function can be re-written in covariant form for 
a general coordinate system as follows: 

dW = gTppdxv, (18) 

with the corresponding wave operator 

r^t = tW* (19) 



associated with the action along a curve, provided 2g^ u = +l v l^ , 
where Y = §^7 a - 

The generalized Dirac equation 

r& - & w 

can be defined along an arbitrary curve and is always covariant. 

Gauge potentials of the form N 1 can be introduced into the system 
by defining 

dx 

P» = = d » w - eA v (21) 

However, in general § A^dx^ ^ and therefore is not an exact differ- 
ential and dsd s W ^ ^jjds but does obey Eqn (5). An analysis of this 
will be given elsewhere. 

As already noted, the above approach deepens our understanding of 
the Principle of Complementarity. The particle properties should be 
directly associated with the path increments. The wave properties 
emerge from Equation (13). 
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• It should be clear that the strict form of the Dirac equation (13) 
pertains to the kinematics and not the dynamics of the motion. It 
describes the kinematics with respect to a local tetrad. The dynamics 
requires further work, which we will do in the next section. Indeed, 
the restriction of the motion to geodesies also explains why we obtain 
distinct energy and momentum levels. Geodesic motion presupposes 
constant momentum and energy, which can either be discrete or form 
a continuum depending on the boundary conditons. Quantum me- 
chanics associates these constants with quantum numbers. 

2.1 Relationship between Quantum and Classical Me- 
chanics 

From a strictly mathematical perspective a particle of rest mass m moving 
in Minkowski space could equally be a cannon ball (a classical object) or 
an electron being ejected by a neutron in beta decay (a quantum mechanics 
object). In effect, the difference between the two objects are determined 
by the initial boundary conditions. In the classical case the line increment 
from which the Hamilton- Jacobi function is derived gives rise to a dual non- 
quantum "wave- function," which is a point mass, given by ip(s) = A5(ms). 
In contrast, in the case of quantum mechanics the same line increment is 
dual to a family of L 2 functions in such way that the initial boundary con- 
ditions coming from the physics are statistical, non-deterministic in nature 
and incorporates quantization. In other words, in the case of a strictly 
classical particle, the mechanics can be determined (in principle) from the 
initial conditions applied directly to the properties of the line increments, 
with the non-quantum "wave-equation" representing a point-mass and not 
contributing any additional information. In the case of a quantum parti- 
cle the opposite appears to be true. It is precisely the "wave-equation" 
that encapsulates the dynamics of the particle, although the solution to the 
"wave- equation" is dependent upon the line increments associated with the 
classical particle. This also gives a new insight into the Principle of Com- 
plementarity. 

Example: Consider a particle of rest mass m moving with uniform velocity 
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u with respect to proper time along the x-axis in Minkowski space. The 
Hamilton- Jacobi function is given by W = m u x x — m it = —m s such that 
Px = fn u x and H = m t. Indeed as a classical particle with x = when 
s = 0, then x = u Q s, where u Q = In terms of the coordinate system (x, t) 
of the laboratory frame this can be written as x = u x t, where u x = {u %) is 
constant. Moreover, this information can also be encapsulated in a family 
of Dirac delta functionals defined by ip{W) = 5 S (W) = W(s) which in term 
of the laboratory frame is equivalent to ip(W(x,t) = 5(x — u x t). In other 
words the wave function as a point mass indicates that at any time t we 
can find the particle in the position x = u x t with probability 1. The wave 
functional per se adds no new information. It essentially encapsulates the 
identical information already obtained by analyzing the Hamilton- Jacobi 
equation derived directly from the line increments ds 2 = u Q dx — tdt. 

On the other hand, when we turn to the problem of decay and an electron 
moving along the x-axis, we have no way of knowing where it will be unless 
we place a detector somewhere enroute, which then measures position at a 
particular instant. Unlike the classical problem, such a measurement does 
not allow us to predict subsequent motion in a deterministic way. Rather, 
we assume that its motion is described by the generalized Dirac equation 
with eigenvector solution ip{W) = Aexp(k(u x — tt)), k a dimension full 
constant. This in turn prompts the question of interpretation. The author 
claims that such interpretations are multiple and depends on the question 
being asked. 

For example, for a detector placed at x Q the arrival times of a beta-decay 
source can be modeled with an exponential distribution #exp(— 9t) with 
mean |. Moreover, for x = x Q , ip Xo (W) = Aexp(2ku x ) exp(—2ktt) defines 
a probability distribution in t. This can be identified with the exponential 
distribution by letting 9 = 2kt = Aexp(2ku x ). In practise, 6 can be 
measured experimentally by considering an ensemble of n beta particles and 
measuring the mean time of detection of each one of them at x . Specifically 
if the ith decay is detecting at time t iy then 

11™ t —t 

^ = i E(^-^-i) = ^ j£ - 

6 nf^ n 
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As a second example, consider n-identical and independent particles in 
Minkowski space with Hamilton- Jacobi function W = J27=i m oSi, (each i 
representing a different particle). Then the function defined by 



iP( Sl ,s 2 ,...,s n ) = n^(si) (22) 

i=i 

= Aexp(kj2(t 2 -Pl-Pl-Pl) (23) 
i=i 

= Aexp(rf)exp ^-fc£(p?+pi + P3)i) (24) 

= Aexp(A;nt 2 )exp ( -kf^Pi J , (25) 



i=l 



where = (p 2 + p 2 , + is a solution to the generalized Dirac equation 
(13), and may be loosely referred to as a "wave-function." As a function in 
Minkowski Space it cannot be normalized, as seen from Equations (23)- (25). 
However, for any fixed t, it can be normalized as an L p function, p > in 
the Euclidean space E 3 . Moreover, as defined it can be used to represent 
either classical or quantum ideal gases according to the initial boundary 
conditions imposed on it. 

p 2 

For example, for a suitable choice of A and k, with = the squared- 
"wave function" can be written as 

^{W) = A(N)^(t)e W (-J-Ylvl^j = A(N)^ 2 (t)e W (-/?£o) 

such that 

^(W( Pl . . . p n ) = A(N) exp ^ p 2 ) = A(N) exp (-fi £ e^j 

is independent of t and defines a Maxwell-Boltzmann distribution for an 
ensemble of n classical particles defined in the center of mass frame. 



16 



In addition, if we also require that the values of pf are such that 

€i = p 2 G {h 2 kf/2m \l = 0, 1, 2,3,...}, E = £ n *h E n i = n > 

and impose indistinguishability conditions (encapsulated in the probability 
amplitude term A) then 

^ 2 (^(Pi • • • Pn)) = M N , ni, n 2, ■ ■ ■ , noo) exp (-p n»e») 

defines a Bose-Einstein statistic or Fermi-Dirac statistic depending on whether 
the occupation numbers n/ range over the set {0, 1, 2, . . .} or the set {0, 1}. 

The above formalism also begs the question as to why L 2 functions are 
needed. Indeed, for the Maxwell-Boltzmann statistics as derived above, it 
would be sufficient to work with any L p function, provided it is also a solu- 
tion to the Generalized Dirac equation. However, when we refer to the Dirac 
equation proper (16) the appearance of % and the imposition of periodic 
boundary conditions to obtain standing waves, requires that the normaliza- 
tion process be restricted to L 2 functions. In reality, quantum mechanics 
is an empirical science with measurements and observations being made in 
real time in the laboratory frame, and any probability interpretation should 
be made with this in mind. In this regard, the above examples highlight 
the importance of boundary conditions when interpreting the significance 
of a wave function tp. In the case of the beta decay problem, the wave 
function ip(x,t) = ip(x)ip(t\x) such that = 1||. In other words, for 

each x we can associate a (conditional) probability distribution with ip{t\x). 
Similarly, in the case of both the Maxwell-Boltzmann and Bose-Einstein 
statistics, ||^(W(pi . . . p„)) = if>(E) = 1||. 



3 Non-geodesic Motion and the Hamiltonian 

In the previous section we related the Hamilton- Jacobi characteristic func- 
tion directly to the solution of the generalized Dirac equation and noted 
that it can represent both classical and quantum solutions of the equation, 
dependent upon the initial conditions. Moreover, from a mathematical per- 
spective if) can be a functional on a space of compact support C£°, or an 
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LP function defined on some domain. However, within the context of Gen- 
eral Relativity, Theorem 1 shows that the general form of the solutions are 
determined only locally and not globally, especially when we consider mo- 
tion along a non-geodesic. Indeed, the existence of non-geodesics suggests 
that other factors other than gravity may be involved. For the purpose of 
quantum mechanics, we will take ip G L 2 (E 3 ), where E is a Euclidean space 
and ^ G L 2 ® H where if is a finite dimensional Hilbert Space associated 
with the spin. For example, ip G L 2 but ip ® £ G L 2 ® H. ^ and ip will also 
be referred to as defining the state of the system, with and without spin 
respectively. 

In this regard, we need also to be aware that from a strict physics per- 
spective, things are more nuanced. Dirac, for example, defines the state as 
the maximum knowledge that we may have about the system. However, 
when one formulates a mathematical theory one is always limited by the 
definitions and restrictions imposed by the "space" within which the theory 
is formulated. For example, in the previous section W was restricted to 
being a Hamilton- Jacobi function which required that momentum be de- 
fined as p a = f^. However, from Lemma 1, we know that ip(kW) is also 

a Hamilton- Jacobi function and in particular = kp a ip(kW), if ip is 

an eigenfunction. Moreover, if k is then chosen to be the complex number 
i, quantum bound states result. Similarly, the Hamiltonian H was defined 
strictly in terms of position and momentum. Consequently, in this context 
the energy states of the system depend only on position and momentum 
and on the imposed boundary conditions. In other branches of physics, the 
state, may also depend on temperature, potential, electric and magnetic 
fields. 

3.1 Hamilton's Equations of Motion 

The Hamilton- Jacobi equation is given by W = J (pdx — Hdt) with the 
understanding that the integration is independent of the path. Using this, 
we now derive Hamilton's equations of motion directly from this definition, 
without making any explicit recourse to the Calculus of Variations. 

First, for each differentiable function W, we write ^ = W. Also note 
that 1 1 ^|| = 1 requires that W = ^ = —H(s), and on taking the deriva- 
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tive of the Hamilton- Jacobi function with respect to s we obtain 

-H(s) = px - Hi, (26) 
with the understanding that for a G {1, 2, 3} 

rir dW TT/ dW dW 

n = ar = ~ H{S) - i^ =p - and "aT = (27) 

Note that these can be written in covariant form = g^f^j- However, for 
simplicity and clarity we will continue to work with local tetrads. 
Differentiating (27) with respect to s gives 

«-*(., and M 

To derive the remaining equation of motion we follow a method introduced 
by Synge and Griffith(see [3]): 

d (dW\ . , 

9 /<9W •. dW \ 

^ + -zrt (30) 



dx a \dx b dt 
d 2 W ■ d 2 W 



x b + — —t. (31) 



Also, 



dx a dx b dx a dt 



d (dW\ . d / ( b dW \\ . 

* = H K'T^M * (32) 



cte a \ dt dx a \ \ ' cte b 



dH.dHJPW. 



cfe a <9p fe dx a dx b 
dH(s) dH d 2 W 



-i. (34) 



^fe Q x aQ x b 

Combining equations (31) and (34) and substituting from (28) yields 

^ ^-« = a (35) 



dx a dx b \ dp b 
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Therefore, 

x b = provided det ( ® W „ - ] + 0. (36) 

dp b \dx a dx b J ^ v ' 

This completes the derivation of the canonical equations of motion. 
In terms of tetrad summation notation these can be rewritten as 

dx a _ h dH{s) df_ h dH{s) 

ds 1 dp b ' ds " 1 dx b • 1 ' 

These are the same equations assumed by Horwitz, Schieve and Piron in 
their work on Stueckelberg theory applied to the Gibb's ensemble[6]. It 
should also be noted, in reference to Eqn. (37) that in tetrad coordinates 
by the principle of equivalence = where D represents the covariant 
derivative. This follows because the affine connection vanishes on a geodesic. 
In terms of a generalized coordinate system both equations of motion can 
be subsumed into the covariant form: 

_ dK Dp* _ dK 

dy> ~~dT ~ ~ 9 d^> (38) 

where = p*u v = p + T* x x >1 x x . Note that in general the covariant 
derivative is necessary on the left hand side, since ^ is not a tensor and 
transforms like acceleration under a change of coordinate systems. On the 
other hand the right hand side of the second equation is already covariant. 



3.2 Non-commutative Operators 

We begin this section by recalling that for this article \1/ = ip l ei represents 
a spinor on a manifold, while ip l are scalar functions (cf. page 2). Also 
if \1> = where £ is a spinor then ip is said to define a scalar field. It 
then follows that d s ip(W) = 7° & ^q^ an d d s ip(H) = 7° 9 ^q^P are operators, 
because of the presence of the 7 matrices, where W is the Hamilton- Jacobi 
function and H is the Hamiltonian of the system. In particular, 

[d 8 ^(W)M(H)]^0. (39) 
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To see this, note that from equations (27) and (28) applied to ifj(W) and 
tp(H) respectively, we obtain 



dil)(W) = 7 1 



a 



dil>(W) 



(40) 



dx a 



and 



dil>(H) 



(41) 



dx' 



Now unless p a = g(t)p a for some function g(t) then they cannot commute. 
Indeed, the condition p a = g(t)p a is equivalent to x a = exp(J* g(u)du) in- 
dependently of a, which defines a non affine parameter along a geodesic 
(cf. [2] Prob., 1.13, n2). Consequently, there do not exist simultaneous 
eigenvectors, and both operators cannot be simultaneously measured, ex- 
cept along a geodesic. This in itself explains some of the difficulty with 
dynamics in quantum mechanics. It is also clear from Equation (5) that the 
non-commutativity is very much related to the presence of the 7 matrices 
and the Dirac Algebra, and although ip(W) or ip(H) may themselves be L p 
functions nevertheless \1/ = £ a spinor, may be interpreted as an element 
of the finite four dimensional Hilbert space associated with the spin of the 
particle. Indeed, the non-commutativity indicates that spin only becomes 
manifest in accelerations, and also suggests a second quantization procedure 
associated with particles moving off geodesies. 

In practice, both equations can be seen as useful depending on the 
circumstances. When the dynamical system corresponds to motion on 
geodesies then p = and has no role to play. In this case the generalized 
Dirac equation (13) can be used. In the event that the dynamical system 
is undergoing accelerations with respect to the laboratory frame then we 
could solve the eigenvector equation 



where = ||/ a || is the norm of the four-force f a . It should also be pointed 
out that since the Principle of Equivalence guarantees the existence of a 
coordinate system x a such that x a = 0, the existence of absolute acceleration 
presupposes the existence of a non-gravitational force such that a a = x 7^ 0. 



dj>(H)Z = 7' 



,0 




(42) 



dx' 
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This also means that in the presence of such an acceleration the rest mass 
will not be a constant but will change according to whether its accelerating 
or decelerating. Accordingly, depending on whether ^ = 0, > 0, < 0, 
three cases arise which could be interpreted as defining a null, time-like, 
and space-like event in momentum space. This, together, with the second 
quantization properties will be investigated in another research paper. 



3.3 Statistical Mechanics and Ideal Gases 

Another interesting application is to theory of Ideal gases. Indeed, as noted 
in Section 2.1 both classical and quantum statistics distributions can be 
derived from solving the generalized Dirac equation, with the difference 
between the two being related to the boundary conditions. In fact, in equa- 
tions (23)-(25) we already assumed the form of the "wave-function" for an 
ideal gas, both classical and quantum. However, it now remains to de- 
rive this fundamental Hamiltonian function as an application of Hamilton's 
equation's of motions. 

First, consider a single particle with a Hamilton- Jacobi function given 
by W = —ms associated with motion along a geodesic in Minkowski space, 
then H = ^ = m. In terms of the equations of motion there exists a 
differentiable function ip = ip(H) such that 

^t- = -ih'p a , where p = (43) 
dx a ds v ' 

which can be re-written in spinor notation as 

= -tW. (44) 



Taking the dot product of (43) with m^- and using the chain rule, (or 
equivalently by taking the inner product \{l a d ^Q^\ l a p a } in (44)) gives 



^ = - WW*. (45) 

To solve, recall that 4>'(H) = d ^}^ and 2p a p a = ^(p a p a ), and consequently 
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H = 2 p a p a - In particular when ip' = kip, an eigenvector, and k an eigen- 
value, solving for ip gives 

rj} = Ae^ Pa . (46) 

If k is real and time dependent then 

ip(t,x) = Aip(t)ip(x\t) = Aexp(kt 2 ) exp(-k(pl + p\ +p 2 3 )), 

is not normalizable. However, for each t, ip(x\t) is an L P (E 3 ) function 
and can be normalized. Indeed, the normalized set of functions {ip(x\t)} 
defines a Markov process in t. Moreover, if A; is a constant as it is for 
geodesic motion, then ip(xi,X2,x 3 \t) = ip(xi, X2, x 3 ). On the other hand, if 
we take k — | then ip(p\t) = ijj(p) is an L 2 function with periodic boundary 
conditions. 

With this in mind, we now extend this result to a system of n indepen- 
dent particles as would occur in some statistical systems, and define the joint 
wave function as an independent product of n single particle eigenfunctions 
with real k, or equivalently as an ensemble of n independent eigenfunctions. 
Written as an L 2 function (although for a classical particle we could equally 
work with L 1 ) 

^ = e ^> ap \ (47) 

Moreover, if the motion is along a geodesic then H = constant and p a p a = m 
for each particle Consequently, 



ip = cexp 



k 

"> -2 -2 ■ 2\ 

-zm^Jt -x 1 -x 2 -x 3 ) 

Z n 



k 



= ei nm . (48) 



Defining T = to be the temperature, where ks is Boltzmann's constant, 
and using the separation of variables for time independent states to write 
ip = ip(t)ip(xi,X2,x 3 ), where 

il>(x 1 ,x 2 ,x 3 \t) = cexp f -£^£ n (x? + x\ + xj) J (49) 
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determines a Maxwell-Boltzmann statistics for free particles at any time 
t. Note also that \ip(xi, X2, xs\t) | 2 defines a normal distribution provided 
the variance is given by a 2 = a 2 = a 2 = a 2 = Equations (47) and 

(48) can be interpreted to mean that the system is in equilibrium, with 
total conserved energy J2 n m - Indeed, for each (local) t the same Maxwell 
distribution occurs. Moreover, if T is not constant then k = k(T) varies 
and in this case, on solving for in (43) one obtains 

Also from equation (44) we should note that ^ is always an exact differen- 
tial, and incorporates k(T). If k = k(T) has no explicit time dependence 
then ip will define a stationary state. Note that in the case of a closed system, 
T is proportional to the population variance of the velocities J2n ^'"^ ■ 



4 Conclusion 

The article has attempted to establish a relationship between the metrics of 
General Relativity and Quantum Mechanics. This has been achieved by first 
relating the metric structure of spacetime to the Hamilton- Jacobi function 
and then using this relationship to derive a Generalized Dirac equation. In 
addition we have derived Hamilton's equations of motion directly from the 
Hamilton- Jacobi equation and then used these equations to determine a dy- 
namical equation for the evolution of the system. The dynamical equations 
derived do not commute with the generalized Dirac equation and conse- 
quently cannot be measured simultaneously. The dynamical equations also 
permit the derivation of the statistical mechanics of the system. Indeed, in 
this paper the Maxwell-Boltzmann distribution was derived directly from 
the equations of motion. 

Finally, it should be noted that we have restricted ourselves to scalar 
fields as defined in the introduction. However, in its most general form, we 
can write 

Tfi? = * (51) 
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where $ would be defined by the physics of the problem. For example, 
Maxwell's equations in Minkowski space can be written in spinor form as 

ia a — = -4tt$, (52) 
dx a K ' 

where <fi = p is charge density, and <p a = j a , a e {1,2,3} is a current 

density. Also in this case, <po = and <p a = H a — iE a , where H a and E a 

are the magnetic and electric fields respectively [7]. Such cases would need 

further study. 
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